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Abstract. The influence of s- and p-wave interactions on trapped degenerate 
one and two-component Fermi gases is investigated. The energy functional of 
a multicomponent Fermi gas is derived within the Thomas-Fermi approximation 
including the s- and p-wave terms of an effective contact interaction. On this basis 
the stability of the dilute gas against mean-field induced collapse due to attractive 
interactions is investigated and explicit stability conditions in terms of Fermi 
momentum and s- and p-wave scattering lengths are derived. Furthermore the 
spatial separation in a two-component Fermi gas is discussed, explicit conditions 
for the onset of component demixing are given, and the density distributions 
in the separated phase are calculated. The findings are summarized in a zero- 
temperature phase diagram for the degenerate two-component Fermi gas. It is 
shown that the p-wave interaction has significant influence on the phase diagram 
of the degenerate Fermi gas and causes new phenomena like absolute stabilization 
against collapse and component separation. It may therefore be useful in the 
context of the envisioned BCS transition in trapped atomic Fermi gases. 



PACS numbers: 03.75.Fi, 34.20.Cf, 32.80.Pj 
1. Introduction 

Looking back at the outstanding experimental achievements in the field of trapped 
ultracold gases composed of bosonic atoms of the last 6 years one wonders what 
exciting new developments one has to envision in the near future One branch 
of these developments will involve dilute and degenerate gases of fermionic atoms. 
Already today several groups ^, |[ ^ were able to cool samples of fermionic atoms, 
^Li or ^"K, to temperatures corresponding to a fraction of the Fermi energy. These 
experiments use either a single fermionic isotope (^"K) in two different hyperfine states 
or a mixture of a fermionic and a bosonic species (^Li/ ""Li) [|[ Q . In this way s-wave 
scatterings between atoms of either the two fermionic components or between bosons 
and fermions allow efficient evaporative cooling of the gas. Remember that in a gas 
of identical fermions s-wave interactions are suppressed by the Pauli principle. 

The experiments on dilute degenerate Fermi gases open a new and unique field 
to study fundamental quantum phenomena. Many of the questions discussed for the 
Bose-Einstein condensed Bose gas also apply to degenerate Fermi gases. Examples 
are the mean-field induced collapse of the dilute gas due to attractive interactions |^ 
or the variety of collective excitations ||] . One of the most challenging experimental 
goals is the realization of a BCS-transition from the normal to a superfluid phase 
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[0, ^. However, a detailed understanding of the "normal" degenerate Fermi gas is a 
prerequisite to achieve that. 

In this paper we investigate the influence of interactions on single-component and 
two-component Fermi gases in the normal phase with special emphasis on the role of 
the p-wave interaction. Therefore in section |^ we construct an effective s- and p-wave 
contact interaction and use it to derive the energy density of a multicomponent Fermi 
gas in Thomas-Fermi approximation. On this basis we investigate the structure of 
one- and two-component Fermi gases at zero temperature. A first question concerns 
is the stability of the dilute gas against collapse induced by attractive interactions. 
In section ^ we derive stability conditions for the one- and two-component gas and 
discuss the influence of the p-wave interaction. A second question deals with the spatial 
separation of the different components in a two-component system. Conditions for the 
onset of component demixing are derived in section |[ We also discuss the structure 
of the density distributions in the separated phase. Finally, in section ^ we summarize 
our flndings by means of a zero-temperature phase diagram for the two-component 
Fermi gas. 



2. Energy functional of a multicomponent Fermi gas 

2.1. Mean-field approximation and effective contact interaction 

In general the description of an interacting quantum mechanical many-body system 
relies on approximations to render the problem solvable. A general procedure to 
obtain different approximation schemes is to restrict the many-body states the system 
can assume to a certain subspace — called model space — of the full Hilbert space. A 
prominent example is the mean-fleld approximation where the many-body state is 
given by a symmetrized or antisymmetrized product of single-particle states. Within 
this mean-fleld model space the many-body problem can be solved; in case of fcrmions 
the solution is obtained by the well known Hartree-Fock equations. 

However, this restriction to a model space of sufficiently simple structure 
encounters a serious drawback. Many realistic interactions between the constituents 
of the many-body system exhibit a strongly repulsive short-range core in addition to 
the attractive part. Typical examples are van der Waals interactions between neutral 
atoms, e.g. the Lennard- Jones potential with an attractive part proportional to r~^ 
and a strong repulsion proportional to r~^^. The repulsive core prevents any of the 
particles to approach each other closer than the core radius, in other words, the two- 
body density distribution is depleted or even zero for particle distances smaller than 
the core radius. These short-range correlations cannot be described within a simple 
model space like the one of the mean-field approximation ^ . 

To resolve this discrepancy we transform the realistic potential into an effective 
interaction suited for the particular model space under consideration. In the case 
of cold and dilute atomic gases the properties of the physical system simplify the 
construction of the effective interaction significantly. This is our starting point: (i) 
We want to describe the properties of a degenerate Fermi gas in the mean-field 
approximation, i.e. the many-body state is an antisymmetrized product state, (ii) 
The system is cold and dilute, i.e. mean particle distance and wavelength of relative 
motion are much larger than the range of the two-body potential, (iii) The system is 
in a non self-bound metastable many-body state, the two-body bound states of the 
potential are not populated. 
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Due to (ii) and (iii) the atoms do not probe the detailed radial dependence of the 
interaction, therefore we can use a simple zero-range effective potential. To describe 
interactions also between particles with relative angular momenta I > the contact 
interaction has to be nonlocal. The mapping of the exact onto the effective interaction 
is done by means of the two-body energy spectrum: We require that the exact two- 
body eigenvalues E^^"^* for positive energies (scattering states) is reproduced by the 
expectation values of the effective interaction calculated with kinetic energy eigenstates 
I nlmi ) 

( I i^h'i' + ''''''' I ) ^ ^nr' (1) 

where q = ^(pi — P2) is the relative momentum of the two particles and m^cd their 
reduced mass. In the following we will consider only particles with equal masses, i.e. 
mred = rn/2. I and m; arc the quantum numbers of relative angular momentum, and 
n is a radial quantum number Guided by (Q) we can construct the operator of 
the effective contact interaction (ECI) explicitly. Up to the quadratic order in the 
relative momentum q, which includes a p-wave contribution, the operator of the ECI 
reads (ITTl 



V 



ECI 



1^ 5(3) (^-) + ^ [^(3) (^-.) . ^2 ^ ^ . ^) ^(3) (^-.) . g.)^ (2) 

mm m 
where r — r/\r\ and h = 1. The properties of the exact two-body interaction enter 
through the s-wave scattering length ao, the s-wave effective volume 69 : and the p- 
wave scattering length ai. They are defined by the first terms of the low- momentum 
expansion of the scattering phase shifts rji (q) of the exact potential 

tan 770 (g) _ 2 tan 7^1 (g) ,3 

Oq — Opg , z oCli- (O 

q (73 J 

The first term of the effective contact interaction is identical to the widely used 
s-wave pseudopotential |l3|. The two additional nonlocal terms are the natural 
extension to describe effective range effects and more important the effects of p-wave 
interactions. 

One should be aware that the effective interaction (||) is designed for the 
application in a mean-field model space. It does not make sense to use this interaction 
in a larger model space or even for the solution of the Schrodinger equation . When 
using erroneously this contact interaction to solve the exact scattering problem one 
does not get back the correct phase shifts. 



2. 2. Energy density of a multicomponent Fermi gas 

We construct the energy density of the inhomogeneous Fermi gas in an external 
trapping potential within the Thomas-Fermi approximation, which was shown [p^ 
to be an excellent approximation for the particle numbers considered here. The basic 
assumption of the Thomas-Fermi approximation is that the energy density of the 
inhomogeneous system is given locally by the energy density of the corresponding 
homogeneous system. 

The calculation of the energy density of the homogeneous multicomponent Fermi 
gas in mean-field approximation is straightforward. The translational invariant part 
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of the Hamiltonian of the trapped gas, consists of the kinetic energy operator and the 
effective contact interaction (|^) 

N N 

^-. = ^E^^?+E-r- (4) 

i i, j>i 

The energy density in mean-field approximation is defined by the expectation value 
of this Hamiltonian calculated with a Slater determinant | ^E* ) 

fhom = -^(*|ifhom|*>. (5) 

The single-particles states | i ) are eigenstates of the momentum operator with 
momenta ki up to the Fermi momentum k. To account for the different mp substates 
that can be trapped simultaneously |^ we characterize the single-particle states by an 
additional component quantum number ^ = 1, . . . , 5 

\i)=\h,^^). (6) 

The many-body state | ) is thus an antisymmetrized product of all single-particle 
states I « ) up to the Fermi momenta of the different components ^. The calculation 
of the expectation value (^) assuming a cubic box of volume V with periodic boundary 
conditions is straightforward. Within the Thomas-Fermi approximation £hom can 
be directly translated into the energy density of the inhomogeneous trapped system 
by replacing the Fermi momenta with local Fermi momenta k^(x) and adding 
the contribution of the trapping potential U^{x). This leads to the general form of 
the energy density of a trapped multicomponent Fermi gas including s- and p-wave 
interactions |l5| 

£[k,,...,ks]{x)^ 4(^^ + ^^^r^J24i^^ + ^^ E 4(^)4(^) 

where the summations range over all components i^, ^' = 1, . . . , S. The interaction is 
characterized by the s-wave scattering length ao, the s-wave effective volume &0j and 
the p-wave scattering length a\. In order to keep the treatment as transparent as 
possible we have assumed equal masses m and interaction parameters for the atoms of 
all components. The local Fermi momentum K^^ix) is connected to the partial density 
p^{x) of the respective component by 

pd^)^-^4('^)- (8) 

For most of the following investigations the detailed form of the trapping potential 
U^{x) does not enter. If, however, the shape of the trap is relevant, then we assume a 
deformed harmonic-oscillator potential 

U{x) = —{ujlxl+Ujjxl+Ujlxl), X= {XI,X2,X3) , (9) 

with a mean oscillator frequency lu — ■^luiuj2'-^3 and the corresponding mean oscillator 
length £ = {muj)~^^'^. 
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2.3. Functional variation and extremum condition 

On the basis of the energy density functional we obtain the density distribution of 
the degenerate Fermi gas by functional minimization of the energy 

E[ki,...,ke]= jd^x E[ki,...,ke\{x). (10) 

with the constraints of given numbers of atoms of the different components 
^=1,...,S 

- jA'x p^{x) = ^ Jd'x 4(x). (11) 

For the variational treatment these constraints are implemented via a set of Lagrange 
multipliers /i^ — the chemical potentials of the different components. This leads to the 
Legendre transformed energy 

F[ni, . . . , ke] = Jd^xJ^[Ki,...,Ks]ix), (12) 

where the associated energy density . . . , ke]{x) is given by 

• ■ ■,Ks]{x) = £[ki, . . .,ke]{x) - ^ (12) 

The set of density distributions p^{x) or local Fermi momenta k^{x) that minimizes 
the transformed energy describes the state of the trapped degenerate Fermi gas. A 
necessary condition for a minimum of the energy F[ki, . . . , kh] is that the derivative 
of the energy density J-[ki, . . . , K'sKx) with respect to the local Fermi momentum 
vanishes for this set 
d 

— — j^J^Iki, . . . , K^](x) =^ for all x, ^. (14) 

OK^(X) 

This condition is also fulfilled by maxima or saddle points, therefore one has to check 
explicitly whether the solutions describe a minimum of F[ki, . . . ,kh]. Inserting the 
energy density (|l^) leads to a set of S coupled polynomial equations for the local 
Fermi momenta at the point x 

+ ^y^^ E [3 4 (^) + 5 4 (^)] ' c = 1, • • • , s. (15) 

These extremum conditions are the starting point for the investigations on the stability 
of the trapped Fermi gas presented in the following sections. 



3. Mean- field collapse 

In the presence of attractive interactions between the trapped atoms, i.e., if the (s- 
or p-wave) scattering length is negative, the dilute gas may undergo a collective 
collapse towards a high-density configuration. The physical origin of the collapse 
is the attractive mean-field generated by the interaction. At sufficiently low densities 
this mean-field attraction is compensated by the positive kinetic energy contribution. 
However, if the density grows the negative mean-field contribution to the energy grows 
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faster than the kinetic energy. At some critical density the kinetic part is not able 
to stabilize the gas against the mean-field attraction and the system collapses. For 
bosonic systems this phenomenon was observed in a series of outstanding experiments 
@. 

The same phenomenon can occur in degenerate Fermi gases with attractive s- or 
p-wave interactions. In the following we derive explicit conditions for the stability of 
one- and two-component systems. 

3.1. Stability of the single- component Fermi gas 

For the single-component system the Legendre transformed energy density ( p^ of the 
multicomponent Fermi gas reduces to 

T,[k]{x) - ^Pix) ~ p\^{x) + 7^^^\x) + (16) 

As a consequence of the Pauli principle the s-wave contact interaction does not 
contribute in a system of identical fermions. However, the p-wave contact interaction 
does contribute and constitutes the leading interaction term in the single-component 
gas. 

To illustrate the origin of the mean-field instability figure ^ shows the energy 
density ( |l6|) in the centre of the trap as function of the Fermi momentum k for a fixed 
value of the chemical potential fi and different p-wave scattering lengths. We express 
all quantities in units of the atom mass m and a length scale £, which later on will be the 
mean oscillator length of the trapping potential. The thin solid curve shows the energy 
density of a noninteracting Fermi gas {ai/i = 0). The minimum defines the actual 
Fermi momentum at the given position x according to the variational description 
discussed above. If we switch on an attractive p-wave interaction {ai/£ < 0), then 
the dominant term in (^) leads to a rapid decline of the energy at large Fermi 
momenta. For sufficiently weak attractions (thick solid curve) a local minimum at low 
Fermi momenta remains. However, if the strength of the p-wave interaction reaches 
a critical value the minimum evolves to a saddle point (dashed curve). Beyond this 
value of the scattering length the energy density drops monotonically and there is no 
low-density solution any more (dotted curve in figure 0) . We conclude that for a given 
chemical potential there is a limit to the scattering length up to which a metastable 
low-density solution exists. 

The analytic criterion for the instability of the Fermi gas is the vanishing of the 
local minimum at low densities; the onset of instability is marked by appearance of a 
the saddle point, where first and second derivative of the transformed energy density 
J-"i(k) vanish. It is advantageous to perform the derivatives of (^6|) with respect to 
density p — K^/(67r^) rather than Fermi momentum. This leads to a condition for the 
existence of the minimum involving the p-wave scattering length ai < and the local 
Fermi momentum ftjf|: 

^ , r .r/-M 2 3(37r)2/3 

~aiK{x)<^— and m[p. ~ U (x)] af < —-^ . (17) 

The second inequality is obtained by inserting the Fermi momentum of the infiection 
point into the extremum condition ( p^ ) to obtain a limit for the chemical potential. 
These are stability conditions for the trapped single-component Fermi gas, in the 

I Alternatively these conditions can be obtained from the extremum condition (|l5|). The collapse 
region is characterized by parameters where has no real positive solution anymore [l5[ . 
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Figure 1. Legendre transformed energy density ^i(k) of the single-component 
Fermi gas as function of the Fermi momentum k for m£^ [fj, — U {x)] = 40. The 
different curves correspond to different p-wave scattering lengths: ai/£ = (thin 
solid line), -0.08 (solid), -0.0915 (dashed), and -0.1 (dotted). 



2.246 



presence of an attractive p-wave interaction ai < 0; any configuration that violates 
these condition will collapse. 

We obtain a more intuitive relation by replacing the Fermi momentum with the 
total number of atoms in a trap of given geometry; here we assume a deformed 
harmonic-oscillator trap (|^) with mean oscillator length £. For a given set of p-wave 
scattering lengths and the corresponding critical chemical potentials according to ( p7| ) 
we solve the extremum condition numerically and obtain the critical particle number. 
The resulting relation between ai and the critical particle number iVcrit of a metastable 
trapped Fermi gas can be parametrized by the stability condition 

(^y) < 1, (18) 

where the prefactor is adjusted according to the pairs of ai and Nd-a obtained 
numerically. 

For the single-component Fermi gas the instability induced by attractive p-wave 
interactions is no severe limitation in the present experimental parameter regime. For 
typical particle numbers of ~ 10^ and a mean oscillator length of £ = 1 /xm the 
p-wave scattering length would have to be of the order of oi — 840aB (where ob is 
the Bohr radius) to cause the instability. However, by exploiting a p-wave Feshbach 
resonance these values of the scattering length seem to be reachable. 

3.2. Stability of the two- component Fermi gas 

The same analysis can be performed for a two-component Fermi gas. In order to 
keep the number of parameters small we assume identical local Fermi momenta 
k{x) — Ki{x) = K2( x) for the two components. A justification of this assumption 
follows in section 11 . For these systems the transformed energy density (|l^) reduces 



to 

JTT^ lUTT^m yTT'^m wn-^m 

where the s-wave effective volume bo is neglected. For this discussion it can be included 
subsequently through a modified p-wave scattering length — s- + bo/ 3. In contrast 
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to the single-component Fermi gas s-wave as well as p-wave interactions contribute. 
Their interplay leads to a rich variety of phenomena. 

As for the single-component system a mean-field instability of the degenerate gas 
can emerge as consequence of an attractive interaction. In the two-component case, 
however, the s-wave or the p-wave interaction may be attractive which causes a more 
complicated stability condition. Evaluating the Fermi momentum associated with the 
saddle point of the energy density ( ]l9| ) leads to the stability condition 

- aoK(f) - 2 [aiK(f)]^ < |. (20) 

It is convenient to discuss the four possible combinations of signs of the s-wave 
scattering length oq and the p-wave scattering length ai separately: 

uq > 0, ai > 0: for purely repulsive interactions the gas is stable against collapse for 
any Fermi momentum. 

ao < 0, oi < 0: for purely attractive interactions both interaction parts cooperate 
and lower the critical Fermi momentum compared to a pure s-wave or p-wave 
attraction. 

ao > 0, fli < 0: the repulsive s-wave interaction stabilizes the system against p-wave 
induced collapse, i.e., the critical Fermi momentum is increased compared to a 
pure p-wave attraction. 

ao < 0, ai > 0: similar to the previous case the p-wave repulsion stabilizes the system, 
i.e. increases the critical Fermi momentum, compared to the case of a pure s-wave 
attraction. Moreover, if the ratio of p-wave and s-wave scattering lengths exceeds 
a limit given by 

^>^«0.311, (21) 
|ao| Stt^''^ 

then the system is absolutely stabilized by the p-wave repulsion. That means, 
despite of the attractive s-wave interaction the mean-field collapse does not occur. 
Notice that the p-wave scattering length required for this absolute stabilization 
is only 1/3 of the s-wave scattering length. 

For repulsive p-wave interactions that are slightly to weak to generate absolute 
stabilization, the density profile of the trapped gas shows a peculiar structure: In 
a region around the minimum of the trapping potential the density is increased 
by typically one order of magnitude compared to peripheral regions. In contrast 
to the outer low-density phase, which is governed by a balance between s-wave 
attraction and kinetic energy, the inner high-density phase is ruled by the balance 
between s-wave attraction and p-wave repulsion. Without the p-wave repulsion 
the high-density phase would collapse immediately. We discuss this phenomenon 
in detail in rcf. pl[ |. 

The last class of interactions (attractive s-wave and repulsive p-wave) is of special 
interest in connection with a BCS phase transition in trapped dilute Fermi gases. The 
formation of Cooper pairs of two fermionic atoms requires an attractive interaction. 
One possibility is an attractive s-wave interaction in a two-component Fermi gas which 
would allow Cooper pairs composed of one atom from each component. The critical 
temperature for the BCS transition increases with the product |ao|K of scattering 
length and Fermi momentum At the same time the mechanical stability of 

the system is a prerequisite. This sets an upper limit to the value of \ao\K — and 
therewith to the BCS transition temperature — described by the stability condition 
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Figure 2. Logarithmic contour plot of tiie critical particle number Af^rit up to 
which the trapped two-component Fermi gas is stable against mean-field collapse 
as function of the ratio of s-wave and p-wave scattering lengths to mean oscillator 
length £ of the parabolic trap. The white area in the upper part indicates the 
parameter region, where no collapse occurs. 



(pO|). However, if an appropriate repulsive p-wave interaction is included, then the 
effect of absolute stabilization can be employed to get rid of this limitation of the 
transition temperature caused by the mean-field instability|§[ 

An overview of the influence of s- and p-wave interactions on the stability of the 
two-component Fermi gas against mean-field induced collapse is given in figure |^. The 
contour plot shows the critical particle number A^crit for a trap of mean oscillator 
length £ as function of the ratios ao/i and ai/£. The particle number is calculated 
numerically in the same way as described for the single-component gas The 
plot confirms the relevance of the p-wave part of the interaction for the properties — 
especially for the stability — of the dilute degenerate Fermi gas. Two major effects 
of the p-wave interaction are: (i) Absolute stabilization of the system by repulsive 
p-wave interactions in the presence of a s-wave attraction (white region) and p-wave 
stabilized high-density phase (region between the dashed line and the onset of absolute 
stabilization), (ii) Mean-field collapse induced by p-wave attraction also for systems 
with repulsive s-wave interactions (lower right quadrant). In addition a competition 
with another effect, the component separation that is discussed in the following section, 
appears (dark gray area). 

As a rule of thumb, the p-wave part becomes relevant as soon as the p-wave 
scattering length exceeds about 1/10 of the value of the s-wave scattering length. A 
detailed discussion of the various effects is presented in ]ll|] . 

4. Component separation 

A second structural transition that can occur in degenerate two-component Fermi 
gases is the spatial demixing of the two components ||l^, |l^ . The physical mechanism 
behind this transition can be understood from the general form of the energy density 

§ Of course it has to be checked that the gain of mechanical stability due to the p-wave repulsion is 
not counteracted by a destructive effect on the Cooper pairs. 
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Figure 3. Contour plots of the local energy density ^2(^1, K2) of a two- 
component Fermi gas in the centre of the trap as function of the Fermi momenta ki 
and K2 for fii = 112 = 250m~^£~^. The left panel shows the case of a moderate 
s-wave repulsion {a^/l = 0.07, ai/l = 0) that still supports identical density 
distributions for both components. The right panel corresponds to a stronger 
s-wave repulsion (ao/l = 0.12, a\/l = 0) which forces the components to demix 
partially. Brighter shading corresponds to smaller values of the energy density. 



of a two-component system resulting from (|T^ 



= ^[('^1(2') -Ml 



1 



207r2m 



)kI{x) 
kHx)] 
kI{x) 4 



- {U2{x) - ^l2) kI{x)\ 
-^4{x)4{x) 

iKf{x)K2{x) 



(22) 



^4{x)4{x)]. 



The driving term for component separation is the s-wave part of the interaction. Its 
contribution to the energy density is proportional to the product of the densities of 
both components pi {x) p2 (a?) oc (x) K2 (^) ■ For strongly repulsive s-wave interactions 
it is energetically favourable to spatially separate the two components such that the 
overlap of the two density distributions and thus the contribution of the s-wave part 
to the total energy is reduced. 

A quantitative example is shown in figure |^. The contour plots depict the local 
energy density T2{ni, H2) for a fixed value of /ii — ^2 as function of the two Fermi 
momenta ni and K2 for two different positive s-wave scattering lengths ao- The left 
panel {uq/I — 0.07) depicts a standard case where the energy density has a single 
minimum at ni — H2 ^ 0. This class of solutions with identical density profiles for 
both components was already discussed in the context of collapse due to attractive 
interactions. 

A new class of solutions appears, if the strength of the s-wave repulsion is increased 
further. As shown in the right panel of figure ^ {ao/£ = 0.12) the minimum at equal 
Fermi momenta has evolved into a saddle point and two energetically degenerate 
minima emerge at ki = 0, K2 > and ki > 0, K2 = 0, respectively. For these 
interaction strengths the density of one component is depleted or even zero and the 
density of the other one is increased in a central region of the trap. 
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4.I. Stability conditions 

The onset of component demixing is indicated by the transition of the minimum at 
Ki — K2 into a saddle point. We can derive a relation that characterizes the onset of 
component separation by looking at the determinant of the Hesse second derivative 
matrix 



D{pi,P2) 



dp{ dpidp2 
9p29pi dpi 



(23) 



where in JF2[ki(pi), ^2(^2)] the Fermi momenta are expressed in terms of the densities. 
Again, it is convenient to use derivatives with respect to the densities, since this 
eliminates the dependence of the transformed energy on the chemical potential 
automatically. Those stationary points of the energy density that are minima (or 
maxima) yield Z3(pi, P2) > whereas saddle points are characterized by D{pi, P2) < 0. 
The parameter sets where the transition from a minimum at p = pi = p2 to a saddle 
point happens are determined by D{p, p) = 0. This is a straightforward generalization 
of the procedure used in section to obtain the stability condition for mean-field 
collapse. Inserting the energy density ( p^ ) into ( |2^ ) leads to 

[2aoK + 4(aiK)^ + 7r][-2aoK + 2(aiK)^ + tt] = 0. (24) 



The first factor resembles the stability condition (20) for mean- field collapse of the 
Fermi gas. The second factor gives a stability condition for component demixing. The 
configuration of overlapping identical density profiles of the two-components is stable 
only if the condition 

ao k{x) ~ [aiK{x)f < ^ (25) 

is fulfilled everywhere in the trap. In regions of the trap where this condition is violated 
the overlapping configuration is unstable against demixing of the components. 

As in the case of collapse it is convenient to consider the different combinations 
of signs of the s- and p-wave scattering lengths separately: 

ao > 0, ai > 0: with increasing s-wave scattering length the critical density for 
component separation decreases. A repulsive p-wave interaction stabilizes, i.e., it 
increases the critical density. If the ratio of p-wave and s-wave scattering length 
exceeds the limit 

^ > ^ - 0-392, (26) 

then the overlapping configuration is absolutely stabilized; no component 
separation occurs at all. 

ao>0,ai<0: the attractive p-wave interaction destabilizes, i.e., reduces the 
maximum density of the overlapping configuration. If the p-wave attraction is 
sufficiently strong such that 

^ > 2(3/7r)2/3« 1.939, (27) 

then the critical density for a mean-field collapse is lower than the critical density 
for component separation. Therefore no component separation but a mean-field 
collapse of the overlapping configuration occurs. 

ao < 0: component separation does not occur, even if the p-wave interaction is 
strongly repulsive. 
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Figure 4. Logarithmic contour plot of tiie critical particle number Af^rit up 
to which the trapped two-component Fermi gas is stable against component 
separation as function of the ratio of s-wave and p-wave scattering lengths to 
mean oscillator length £ of the parabolic trap. The white area in the upper part 
indicates the parameter region, where no separation occurs. 



The twisted behaviour of the p-wave interaction with respect to component separation 
can be understood from the structure of the energy density (|2^). The driving force 
of component separation — the s-wave interaction — contributes a term proportional to 
the product of the densities of both components. Thus for repulsive s-wave interactions 
a configuration with fully overlapping density distributions is higher in energy than 
a separated configuration with small overlap but higher absolute densities. The p- 
wave interaction, however, has a term [pi p^l^ -I- fh/^ p-;^ proportional to the product 
of the densities to some power and an additional term that depends on the sum of the 
densities, [p^^"^ + ^2^"^]- Thus a strongly repulsive p-wave interaction disfavours both, 
large overlaps of the two density distributions and large densities of a single component. 
Therefore a p-wave repulsion alone cannot induce component separation. Moreover, 
it stabilizes against separation in connection with a repulsive s-wave interaction. 

A pictorial overview of the influence of s- and p-wave interactions on the onset 
of component separation is given by the stability map in figure ^. The contour plot 
shows the logarithm of the critical particle number iVcrit as function of the s- and p- 
wave scattering lengths. This critical particle number is determined from the solution 
of the extremum condition (Ha) that reaches the critical Fermi momentum (Ga) in 



the trap centre. Like in section 3.2 a deformed parabolic trap of mean oscillator 
length i is assumed. As soon as the particle number N = Ni = N2 exceeds the 
critical particle number A^crit a demixing of the components appears in the trap centre. 
The plot clearly reveals the influence of the p-wave interaction: Attractive p-wave 
interactions, on the one hand, reduce the critical particle number. Moreover, for those 
interaction parameters in the dark gray area the p-wave attraction is so strong that 
the mean-field collapse of the overlapping configuration happens before the critical 
density for component separation is reached. Repulsive p-wave interactions, on the 
other hand, stabilize the overlapping configuration against component separation and 
prevent separation completely for interaction parameters in the white region in figure 

i 
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4-2. Structure of the separated configurations 

With a few examples we want to discuss typical shapes of the demixed density 
distributions that appear beyond the critical particle number. To discuss the general 
phenomenology we restrict ourselves to pure s-wave interactions [a\l i — 0) and equal 
particle numbers N = Ni = N2 in a spherical symmetric trapping potential. 

Even for fixed particle numbers and interaction strengths there exists a rich 
variety of different density distributions that may appear in case of component 
separation. In general it depends on the detailed dynamics of the system which of 
the energetically (almost) degenerate states is realized. We can distinguish two major 
classes of separated solutions: Those with equal chemical potentials /ii = /i2 for both 
components and those with /^i /^2■ 
For equal chemical potentials III = iJ,2 [more precisely Hi — Ui{x) = /^2 — ^^2(2?)] the 
extremum condition ([l^) is symmetric with respect to exchange of the two components. 
If the s-wave scattering length is large enough for phase demixing then there are two 
energetically degenerate solutions {p^'^-'(x), P2^'(x)} and {p\^\x), P2^\x)} (compare 
figure where 

p[''\x) = )(f), pi^\x) = p^^\x). (28) 

In general to construct density profiles with equal number of particles in both 
components we can switch between the two solutions depending on the position in 
the trap. The most simple ansatz is to take solution (A) in one hemisphere and (B) 
in the other 



\yr^ + z^) for z < 

which guarantees A^i = A^2- Examples for the density profiles p^^'^^ir^ z) along the z- 
axis obtained for different s-wave scattering lengths are shown in figure |. The upper 
left panel corresponds to an interaction strength that is slightly to weak to cause 
separation. A tiny increase of ao j I leads to the onset of component separation (upper 
right panel), a small density wiggle appears on top of the smooth density profile in 
the trap centre. Further increase of the s-wave repulsion enhances this demixing and 
finally leads to a nearly complete spatial separation where one component is localized 
in one hemisphere and the second one in the other hemisphere (lower right panel). 

The ansatz (^9|) is by no means unique. One could subdivide the trap volume 
into several domains and switch between the two solutions in each of them. For 
example, in a spherical trap one could imagine domains that look like pieces of a 
cutted pie. Within the Thomas-Fermi approximation the total energy is not affected. 
However, if we go beyond the Thomas-Fermi approximation, e.g. by including gradient 
corrections to the kinetic energy, the strong density variations at the interface of the 
domains produce a positive energy contribution. For typical particle numbers of 10^ 
the relative contribution of the gradient corrections for a cutted distribution like ( p9|) 
is in the order of 10^^ Even though this is negligible, we assume that the 

configuration with minimal interface is preferred — the axially symmetric profile ( pof ) 
in this case. 

There is also a class of spherically symmetric solutions with A'^i — N2 that was 
discussed by several authors |l|, |l^. However, it turns out that these spherical 
solutions are higher in energy than the axial ones. In figure ^ two examples for 
these spherical configurations are shown. The plots also contain the respective values 
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Figure 5. Axially symmetric configuration: density profile = 0,z) 

along the z-axis for two-component Fermi gases with A^i = JV2 = 10^ and 
different s-wave scattering lengths ao /£■ The distribution of the second component 
paxial^y. _ Q^^^ _ ^axial^j. _ 0,-2) is obtained by inversion of the ordinate. An 
animated GIF of the density distributions is available from the article's abstract 
page in the online journal fed/. 
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Figure 6. Spherical configuration: density profiles p^^^ir = 0,z) (solid line) and 
p'^^{r = 0,2) (dashed line) along the 2-axis with A''i = N2 = 10^ for different 
s-wave scattering lengths ao/l. 



of the total energy in units of the energy Eq — -s/GN of the noninteracting two- 
component gas with the same particle number. For comparison the energies of the 
axially symmetric configurations are given in figure |^. It should be noted that in the 
spherical case the components have slightly different chemical potentials. 

In an equilibrium situation or during a quasi-stationary process that leads from 
overlapping density distributions pi{x) = P2{x) to a separated configuration, e.g. by 
slowly increasing the s-wave scattering length using a Feshbach resonance ||l^, we 
expect the energetically lowest, i.e. the axial configuration with least interface shown 
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Figure 7. Zero-temperature phase diagrams of a two-component Fermi gas in the 
parameter plane spanned by s-wave scattering length ao/£ and particle number 
N = Ni = N2. The three diagrams correspond to different values of the p- 
wave scattering length: (a) ai/i = 0, (b) ai/i = 0.02, (c) ai/i = -0.02. The 
white area marks the parameter region where both components show identical 
overlapping density profiles. The dark gray area indicates the region of mean-field 
collapse and the light gray region characterizes parameter combinations where 
both components separate (see text). An animated GIF of the phase diagram is 
available from the article's abstract page in the online journal fef]/. 



in figure ^, to be realized. This evolution of the density distributions is illustrated by 
an animation that is available from the article's abstract page in the online journal or 
from 1^. 

5. Summary: Phase diagram of trapped degenerate Fermi gases 

To summarize our findings on the stability against mean-field collapse and component 
separation we construct a phase diagram of the two-component Fermi gas at zero 
temperature. The relevant parameters are the s- and p-wave scattering lengths 
and the particle number N = Ni — Figure shows the phase diagram in 

the plane spanned by the s-wave scattering length ao/i and the logarithm of the 
particle number N for three different values of the p-wave scattering length ai/£. 
We can identify three classes of density distributions — or phases — of the trapped two- 
component Fermi gas: identical overlapping density distributions for both components 
(white region), partially or completely demixed density profiles (gray region), and the 
p-wave stabilized high-density phase (light gray area) . In addition the collapse of the 
system due to an attractive mean- field can occur (dark gray). 

In absence of the p-wave interaction the phase diagram has a rather simple 
structure shown in part (a) of figure ^ If the s-wave interaction is repulsive (oq > 0) 
the system enters the separated phase if the quantity \/Nao/i exceeds a limiting 
value. Similarly the gas collapses in the presence of attractive s-wave interactions if 
the quantity -v^iV oq /(■ exceeds a another maximum value. 

The inclusion of p-wave interactions modifies this picture significantly. Panel (b) 
of figure ^ depicts the phase diagram for a rather weak repulsive p-wave interaction 
with ai/i — 0.02. For s-wave interaction strengths below a certain limit the 
overlapping phase extents up to very high particle numbers; the p-wave repulsion 
generates absolute stabilization against collapse and component separation. Moreover, 
the p-wave stabilized high density phase emerges near the threshold for absolute 
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stabilization against collapse (light gray region). 

Dramatic effects on the phase diagram are caused also by weak attractive p- 
wave interactions as shown in part (c). The maximum particle numbers possible in 
the overlapping phase are significantly reduced, i.e. the p-wave attraction promotes 
collapse as well as component separation. In addition it causes the collapse of the 
separated configuration, an effect that is in close relation to the collapse of a single 
component Fermi gas. If the density or Fermi momentum of one of the components 
in the separated configuration exceeds the critical value given by (|l^) then this 
component will collapse. Actually this causes a distinct reduction of the parameter 
region where a stable separated configuration occurs. 

Our investigations show that the p-wave interaction can have big influence on the 
properties of trapped degenerate Fermi gases. In the single-component system it is the 
origin of the mean-field instability. In a two-component systems it modifies the phase 
diagram substantially and gives rise to completely new phenomena like the absolute 
stabilization against s-wave induced collapse and separation. In this way the p-wave 
interaction may be very useful for the envisioned transition to a superfluid state in 
trapped dilute Fermi gases. 
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